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Abstract 

We show that the nonlinear massive gravity model of de Rham, Gabadadze, and Tolley admits 
exact plane gravitational wave solution whose waveform obeys the two-dimensional Helmholtz 
equation. The solution is valid for arbitrary values of the coefficients of the cubic and quartic 
terms. In the massless limit the solution reduces to the Aichelburg-Sexl metric in general relativity. 
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I. INTRODUCTION 



the Fierz-Pauli model of massive gravity 



Massive gravity theories have been of interest for the past decades. Various reasons for 
this interest range from pure field theoretic considerations to the recently raised problem 
of the accelerated expansion of the universe. In fact it has been recently suggested, e.g. 
in Refs. [l, 2|, that massive gravitons may be considered as an alternative to dark energy, 
explaining the accelerated expansion. 

The pioneering work by Fierz and Pauli jsj resulted to a now well-known model of mas- 
sive gravity. That model with its unique property of being the sole massive extension of 
the general theory of relativity which is linear, is faced with subtle theoretical problems 
such as the extra helicity degrees of freedom and the so-called van Dam-Veltman-Zakharov 
discontinuity (4, 5] which implies that the theory does not lead to the general relativity in 
the zero- mass limit. It has been attempted to remedy this by adding nonlinear terms to 

6|. However they failed to do so as shown by the 
analysis of Boulware and Deser [7| according to which, an unwanted peculiarity is always 
associated with these models: the appearance of ghosts and its consequent instability. 

A great deal of progress has been recently made in this field by the new massive gravity 
(NMG) model proposed by Bergshoeff et.al. which have proved to be unitary and ghost 
free jol, [loj], and there are claims that it is renormalizable However, the latter model 
is restricted to three space-time dimensions. A recent review of the theoretical aspects of 
various massive gravity models may be found in Ref. l^j], see also Ref. (l^ . 

Recent works on massive gravity in the context of effective field theories [lj], [l^ have 
resulted in the hope that there might be some way out of this dilemma. In a recent paper by 
de Rham et. al. 16| a four- dimensional nonlinear massive gravity model (the dRGT model) 
was proposed which is free from the Boulware-Deser ghost up to and including quartic order 
away from the decoupling limit. The dRGT model is essentially constructed by algebraic 
resummation of the nonlinear terms in a previously proposed effective field theory of massive 



gravity 
in Ref. 



7H19j]. Some related works may be found in Refs. 2CH23| . The model was used 



24j to explain the cosmic acceleration from 0- helicity gravitons. The Vainstein 



mechanism of the model was investigated in Ref. |25lj. The effect of massive gravitons on 
the stell ar g ravitational field was studied in Ref. [26|. Strong coupling effects was studied 
in Ref. 27| . Exact Schwarzschild-de Sitter black hole solutions of the dRGT model have 
been presented in Ref. (28 1 and they have shown there to go smoothly to the Schwarzschild 



and Reissner-Nordstrom solutions of general relativity in the zero-mass limit for a specific 
choice of the parameters. 

In the present work we seek exact plane gravitational wave solutions in the context of the 
dRGT model. This is well motivated by the current interest in gravitational waves in general 
and in exact plane wave solutions, in particular. For example, it is believed that certain 
classes of plane gravitational waves, namely gravitational shock waves, have dominant role in 



scattering process at ultrahigh-energy regime 



29 



30J . Also, the propagation of gravitational 



waves from various sources have been studied in order to place an upper bound on graviton 
mass, see e.g. Ref. 3jJ and references therein. 

In the next sections, we shall start with a brief review of the dRGT model and by solving 
the equations of motion, we show that the model admits exact plane gravitational waves. 



II. THE DRGT MODEL 

Let us represent the Minkowski metric in a chart {x a }, a = 0, 1, 2, 3 by 

7/» = Vabd^d^ (1) 

where r\ ah = diag(—l, 1, 1, 1), and (f> a = x a are the Stiickelberg fields in the unitary gauge. 
The space-time metric is then incorporated via 

i\ = g» x ixv (2) 



The dRGT model of massive gravity is described by the following Lagrangian density 
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£ = lM* v ^g-(R-m 2 U) (3) 

where the first term in the parentheses is the curvature scalar which together with the 
coefficients outside constitute the Einstein-Hilbert Lagrangian density. The second term 
defines the mass term and is given by 

U = a 2 C {2) + a 3 £ (3) + a 4 £ (4) (4) 

with a 2 ,a 3 , and being dimensionless coupling constants. The sign of a 2 is negative, 
otherwise tachyons are allowed in the model. In fact, one may set a 2 = — 1 without loss of 
generality. However, in what follows we also consider the effect of a 2 having the (wrong) 
positive sign. 



The explicit form of £( 2 > 3 ' 4 ) can be expressed in terms of the tensor kP v defined by 

k\ = Sj£- P" v (5) 
where (3 fl u satisfies P^x^u = l^v Defining 

(/C n ) = g/j,i/K^ ai K 1 a2 ■ ■ ■ K n 1 (6) 

we have 

£( 2 ) = (.) 2 -(. 2 ), (7) 

£( 3 ) = («) 3 -3(k}(« 2 }+2(« 3 }, (8) 
£( 4 ) = ( k ) 4 _6(/c) 2 ( K 2 )+3(k 2 ) 2 

+8(k)(k 3 )-6(k 4 ). (9) 

The equations of motion associated with the above Lagrangian read 

v = _J* J(*yr=9W) (10) 
v-g 09 1» 

where i = 2,3,4. On the other hand, from Eqs. ([Z])-© we have the following relations 

5£W 



kT - H 1 " + (n)^ - k^), (11) 

3 - 3 
5 9llv 2 2 

5C^ 



8g t 

<5£ (3) 3 3 

-K» a H av + -H" a K av + 3H" U - 6kT 



+ -C®(gV - O - 3{k)(H^ - (12) 
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&9, 



hp 



12H» a H av - 2AH tiV + 48^ 
-18^ a H au - 18H^ a K au 
+2£ (3) (c^ - k» v ) - 6£ {2) {H^ - & 
+2(k) {3n\H au + 3H» a K au 

+6H^ U - 12kT) (13) 

where H^ u — 8£ — 7^ v . By inserting these into Eq. ffTUl) . the explicit form of the equations 
of motion are obtained. The resulting expressions are algebraically complicated and not 
illuminating. In the next section we present the relevant relations for the case of a plane 
gravitational wave. 



III. GRAVITATIONAL WAVE SOLUTION 



We consider a plane gravitational wave ansatz given by the following line element 

ds 2 = —dudv — F(u, x, y)du 2 + dx 2 + dy 2 (14) 

in which u = t — z and v = t + z are the null coordinates and F(u, x, y) is to be obtained from 
the equations of motion. By taking F(u, x, y) equal to the product of arbitrary functions of 
u and solutions of the Laplace equation (in terms of x and y), this is an exact gravitational 
wave solution in general relativity. In this chart the Minkowski metric is given by 

/ — | \ 



and hence, 



-\ 

10 

\ o o o i / 

( 1 o o o\ 

-2F 1 

10 

\ 1/ 



(15) 



(16) 



from which we obtain 



(17) 



/ 1 o o o\ 

—F 10 
10 
\ 1/ 

There are several other square roots, but they are incompatible with the above ansatz. We 
thus have 



This results in 



/ o\ 
F 

\0 0/ 

£ w = 0, 



(18) 



6£W 
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= 0, 
= 0, 



(19) 
(20) 
(21) 



8tf®_ 



and 

/o o\ 

2F 
OOOO 
\0 0/ 

Inserting these into the field equation (fTOj) and noting that for the metric given in Eq. (H 
we have 

/o o\ 

-2V 2 F(u,x,y) 


\0 0/ 



(22) 



(23) 



where V is the Laplacian in the transverse plane, we obtain 



V 2 F(u,x,y) = a 2 m 2 F(u, x, y). 



(24) 



Thus we obtain a plane wave solution provided the waveform satisfies a two-dimensional 
Helmholtz or modified Helmholtz equation, depending on the sign of a 2 . Therefore we 
obtain 

F{u,x,y) = frWFtfay) + f 2 {u)F 2 {x,y) (25) 

in which fi(u) and f2{u) are arbitrary functions of u and Fi(x, y) and F 2 (x, y) are solutions to 
the two dimensional (modified) Helmholtz equation. This describes a plane wave propagating 
in the ^-direction with the speed of light. The solutions to the (modified) Helmholtz equation 
can be expressed, in polar coordinates, in terms of the (modified) Bessel functions. The 
singularity at x = y = then suggests a source located there. Thus the above plane wave 
corresponds to the gravitational wave generated by that source. In particular, this may be 
compared with the Aichelburg-Sexl solution js^] 



F(u, x, y) = C5(u) ln(x 2 + y 2 



(26) 



with C being a constant proportional to the particle momentum. This describes the general 
relativistic shock wave due to a massless point-like particle propagating longitudinally. 

The modified Helmholtz equation is usually associated with diffusion-like phenomena and 
hence the non-negative values of a 2 are also disregarded from this point of view. 

One can study the physical properties of the plane wave solution Eqs. f|T4l and fT25l) . 
by considering its effect on the motion of test particles. This can be achieved by solving 
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the geodesies of the space-time for the trajectory of test particles or by taking the geodesic 
deviation equations into account. This has been considered for the case of plane waves 



of type given by Eq. (1261) in several works including Refs. [3314351] . Since the equation for 
geodesies involves the first derivatives of F(u, x, y) with respect to the transverse coordinates 
(through the connection coefficients), the effect of the graviton mass, the right-hand side of 
Eq. ( 124)) . can then be traced from the terms containing such derivatives. 

Now it is interesting to ask if the above plane wave, Eq. (|14p . is a solution to massive 
gravities with a general nonlinear potentials given by 

U gen = -C^ + Y, C fen (27) 



i=3 



in which 



4?n = (H) 3 + Cf (H) (H 2 ) + Cf (tf 3 > , (28) 
4i = C?\H)* + C^(H) 2 (H 2 ) + C^(H 2 ) 2 

+ci 4) (H)(H 3 )+C^(H±), (29) 
4ei = Cf\Hf + C^(H) 2 (H 3 ) + ci 5 \Hf(H 2 ) 

+Ci 5) (H)(H 4 ) + Ci 5) (H 2 )(H 3 ) 

+C<i ] (H){H 2 ) 2 + cf ) (# 5 ) (30) 



and similar expressions for higher order terms. Here (H n ) is the trace of (H^ u ) n , and C 



(n) 



3 



are generic coefficients. Now by using the ansatz given by Eq. (1141) . we can repeat the 
previous calculations to obtain 

41 = 0, (31) 



-gen 



0, (32) 



Thus the plane wave is a solution to Eq. ( flOl) for massive gravities with a nonlinear potential 
with non-tuned coefficients. However, such models are not free of ghosts in general. 



IV. CONCLUSIONS 



We have shown that the dRGT model of massive gravity admits plane gravitational wave 
solutions. The waveform only depends on the quadratic term in the interaction potential and 
obeys the Helmholtz equation for negative values of 02. It reduces to the general relativistic 

7 



result, the Laplace equation, for a vanishing 02. Thus we obtained the generalization of the 
Aichelburg-Sexl solution in the framework of the dRGT massive gravity model. The graviton 
mass enters the solution via the massive Helmholtz equation. The solution is valid regardless 
of values of the coefficients in front of the cubic and quartic terms in the lagrangian. In the 
massless limit the solution is the same as in general relativity. The plane wave solution may 
be used to study the scattering of particles in the gravitational field of a massless particle 
in the context of massive theories of gravity. 
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